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We describe a controllable and precise laser tweezers for Bose-Einstein condensates of ultracold 
atomic gases. In our configuration, a laser beam is used to locally modify the sign of the scattering 
length in the vicinity of a trapped BEC. The induced attractive interactions between atoms allow to 
extract and transport a controllable number of atoms. We analyze, through numerical simulations, 
the number of emitted atoms as a function of the width and intensity of the outcoupling beam. We 
also study different configurations of our system, as the use of moving beams. The main advantage 
of using the control laser beam to modify the nonlinear interactions in comparison to the usual way 
of inducing optical forces, i.e. through linear trapping potentials, is to improve the controllability 
of the outcoupled solitary wave-packet, which opens new possibilities for engineering macroscopic 
quantum states. 



PACS numbers: 42.65.Jx, 42.65.Tg 



I. INTRODUCTION 



The achievement of Bose-Einstein condensation (BEC) 
in dilute gases of alkali atoms 1] has driven the re- 
search on the design of new tools for the manipulation 
and coherent control of atomic ensembles. In the last 
years, an intensive study of different mechanisms for this 
purpose has been carried out both theoretically and ex- 
perimentally. Among the most important contributions 
we must cite the realization of atom mirrors guides 
[1, the design of atomic accelerators both in linear 
and circular geometries [1, @, 0] ■ Atom lasers have also 
been developed, first based on the use of short radio- 
frequency pulses as an outcoupling mechanism, flipping 
the spins of some of the atoms to release them from the 
trap [3]. Later, other coherent atomic sources were built 
leadin g to pulsed, semicontinuous or single-atom lasers 

The control of coherent atomic beams is a challenging 
problem in physics due to its potential applications in 
multiple fields like atom interferometrv | 15i| . superposi- 
tion of quantum states [l^ , atom clocks [l7| |. or quantum 
information fl^, among others. Many of these devices 
take advantage of nonlinear interactions between atoms, 
which are ruled by the value of the scattering length 
a. The adequate tunning of this parameter has been 
possible with the use of Feshbach resonances [l^ and 
has yield impressive effects like the macroscopic collapse 
of matter waves [lO] or the creation of atomic solitons 
[21I . [2^ . [23j . The recent realization of optical control of 
Feshbach resonances ^23] has paved the way for the exper- 
imental demonstration of many theoretical proposals on 
nonlinear waves in Bose-Einstein condensates with spa- 
tially inhomogcneous interactions [25', '26*1 including the 
dynamics of solitons when the modulation of the nonlin- 
earity is a random ^ , linear [23| , periodic [11] , localized 
function (s^l or step-like function [sij . 



In this paper, we propose the use of spatially- 
dependent scattering length as a tool for designing pre- 
cise atom tweezers which are able to extract a portion of 
atoms from a BEC. Other methods have been proposed 
for this purpose, as in Ref. [32] and Ref. [33] in combina- 
tion with spatial light modulators for splitting the atomic 
cloud. Our device is inspired on a coherent atomic source 
based on the spatial modification of the scattering length 
[3^ . [3^ . which produces a highly regular and control- 
lable number of atomic pulses by modulating a along the 
trapping axis of a BEC. As we will show in this paper, 
in comparison with linear traps which do not alter the 
value of a, spatially-tailored nonlinear interactions yield 
robust control of the number of atoms that can be ex- 
tracted from a BEC reservoir, providing new ways to the 
creation of macroscopic superposition of quantum states 
from arrays of Bose-Einstein condensates [1^ . 

The structure of this paper is as follows: in Section 
im we introduce the configuration of the system and the 
mean field theory which is used in this work. In Section 
lllll we study by means of numerical simulations the mech- 
anism of emission by employing a static laser tweezers. 
We also analyze the number of atoms emitted as a func- 
tion of the main parameters of the system. Finally, in 
Section HVl we analyze the use of moving tweezers which 
trap atoms by crossing the condensate, and compare our 
results with the linear case in which there is no spatial 
variation of the scattering length. 



II. SYSTEM STUDIED AND MODEL 
EQUATIONS 

We will assume that a BEC is strongly trapped in the 
transverse directions (a;, y) and weakly confined along the 
longitudinal one (z) leading to a cigar shaped configura- 
tion. We will consider the effect of a spatial variation of 
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FIG. 1: [Color online]. Sketch of the system we will study 
in this paper for the case of optical tweezers which are used 
to extract a given number of atoms from a BEC reservoir. 
A rotating mirror is used to move a laser beam over a BEC 
which is transversely trapped by magnetic confinement, and 
optically along the direction of motion of the laser. 



the scattering length along z which can be switched from 
positive to negative values by the optical control of Fes- 
hbach resonances by means of a laser beam. The region 
of negative scattering length can be varied in size and 
displaced along z by simply focusing or moving the laser 
beam. The choice of an optical control ^24| instead of a 
magnetic one [s^l allows for a faster and easier manipu- 
lation of the spatial variations of the scattering length. 

The mean field description of the dynamics of the BEC 
is provided by a Gross-Pitaevskii (GP) equation of the 
form: 
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A2^' + F(f)1' + C/|1'p*, 



(1) 



where is the order parameter, normalized to the num- 
ber N of atoms in the cloud: TV = / l^pd^r. V = 
A-K?i?a/m characterizes the 2-body interactions deter- 
mined by the value of the scattering length a. The cloud 
of N equal bosons of mass m is tightly trapped in (a;, y) 
by a harmonic potential Vj^ of frequency and weakly 
confined along z by the effect of an optical dipole trap Vz 
that can be produced by a laser beam of a given width 
along z [H, l39j . The mathematical expression for the 
potential is thus: 
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where Vd is the depth of the (shallow) optical dipole po- 
tential and L its characteristic width along z. To fix 
ideas, we will present specific numbers in this paper cor- 
responding to ^Li, using the experimental parameters of 



Ref. (24I, where Vd ~ hi'±, being iy± = 1 kHz the fre- 
quency of the trap in the transverse plane {x,y), which 
yields to a transverse radius r± w 3fim. The other nu- 
merical values used in our simulations are L = 4rx, 
= 3 • 10^, Wc = 5Ar± and a = -1.4 nm. We must 
stress that our basic ideas should hold for different atomic 
species like ®^Rb and ^"^^Cs with an adequate change of 
the parameters used. 



III. ATOM EXTRACTION WITH STATIC 
TWEEZERS AND CONTROL OF THE ATOMIC 
WAVEPACKET 

The problem we will face is the controllable extraction 
of atoms from a BEC reservoir. We will consider a system 
configuration in which a trapped cigar-type condensate is 
partially overlapped by a laser beam in a similar configu- 
ration to the one described in Ref. [s^l ■ Under adequate 
conditions, the laser changes the local value of the scat- 
tering length in part of the cloud. If a is locally switched 
to large enough negative values, a burst of atomic soli- 
tons can be emitted from the condensate. If the distance 
between the laser beam and the center of the cloud is 
kept constant, a emitted soliton may rebound inside the 
laser beam and thus remain trapped out of the reser- 
voir [sB]. Once the atoms are extracted, the laser can 
be moved away from the condensate. This also allows to 
control the position of the soliton in the z direction. This 
idea of extracting atoms is radically different than using a 
usual dipole trap to extract atoms without switching the 
scattering length to negative values. In the later case, 
the atoms will perform Josephson oscillations between 
the reservoir and the tweezers and it is only possible to 
extract a significant portion only at times exactly match- 
ing the maxima of the periodic motion. Thus, the role 
of nonlinear interactions is essential in this static config- 
uration to guarantee that once extracted the atoms will 
not go back to the reservoir. 

Our results are based on numerical simulations of Eq. 
([T]). All results to be presented in what follows have 
been obtained using second order in time split-step pseu- 
dospectral solvers, with the spatial derivatives being eval- 
uated by Fourier methods [4C|. In Figs. [2{a) and[D(b) we 
show some numerical simulations showing how our non- 
linear tweezers work. In both pictures, we have plotted 
pseudocolor images of the cloud density. The horizontal 
axis is z and vertical axis is time. The extraction of atoms 
is made with a gaussian-profile laser beam. If the scat- 
tering length is switched to large-enough negative values, 
part of the cloud is extracted. Once the atoms have left 
the reservoir, the laser beam is moved away from the con- 
densate dragging part of the atoms. In Fig. ^a.) the laser 
beam extracts one soliton and controls its position along 
z. In Fig. Hfb) two solitons are extracted with two differ- 
ent laser beams and their paths joined at a given point. 
Time in the vertical axis goes from t = to t — 500 /i^± 
in Fig. H (a) and from t = to t = 1800/j/± in Fig. 
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FIG. 2: [Color online] Controlled emission of atomic solitons 
from a BEC reservoir. The extraction was made by employing 
a gaussian-profile laser beam. Once the soliton is emitted, the 
beam is separated from the condensate dragging the emitted 
atoms. In (a) one laser beam extracts one soliton and controls 
its position in z. In (b) two solitons are extracted with two 
different beams and their paths joined. The range of times, 
shown in the vertical axis is from t — to t = 500/:^^ in (a) 
and to t — 1800/!^x in (b), where i'± — 1 kHz is the radial 
trapping frequency. The horizontal axis is 60 times the width 
L of the optical dipole trap that confines the condensate in the 
z direction. The figures at the right display the profiles of the 
condensates and emitted solitons (in black continuous lines), 
and of the laser beams (in dotted black lines and shaded) for 
different times. 



Wljo) J where v± = 1 kHz is the radial trapping frequency. 
The horizontal axis spans 60 times the width L of the 
optical dipole trap that confines the condensate in the 
z direction. The figures on the right represent the atom 
density showing the profiles of the reservoir and the emit- 
ted solitons (in black continuous lines), and of the laser 
beams (in dotted black lines and shaded) for three differ- 
ent times. As it can be appreciated our control method 
allows a robust control on the extracted atoms. 

Let us now consider a slightly different configuration 
in order to show the robustness of the method. We now 
address the case of a laser beam which is more intense 
and narrow than in the previous simulations. In this case 
a high nonlinear interaction stripe is generated and some 
atoms will be attracted to this thin region of negative 
scattering length and will be trapped. As in the previous 
case, the position of the extracted atoms along z can be 
controlled by moving the laser. In Fig. [3)^ a) and Fig. [3)^b) 
we show a similar representation as in Fig. [Ifa) and Fig. 
[2ljb). In this case, a much narrower beam has been used 
in order to suppress internal rebounds of the extracted 
atomic beam. In Fig. [3)Jb) two matter waves are emitted 
employing two different laser beams. By displacing the 
lasers, it is possible to control the relative position of the 
extracted solitons. Vertical axis corresponds to times in 
the range from t = Otot = 5000/j^x for both figures. The 



plots at the right display the profiles of the condensate 
and emitted solitons (in black continuous lines) and of 
the laser beams (in dotted black lines and shaded) for 
three different times. 
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FIG. 3; [Color online] Same as in Fig[T] for a narrower and 
more intense laser beam. Vertical axis corresponds to time 
spanning the interval from t = to t = SOOO/z/^. The other 
parameters are as in Fig. [2] The figures at the right display 
the profiles of the condensates and emitted solitons (in black 
continuous lines), and of the laser beams (in dotted black lines 
and shaded) at three different times of propagation. 



The number of extracted atoms depends on the laser 
intensity. In Fig. 0] we can see the percentage of atoms 
attracted and trapped by the laser tweezers as a function 
of the intensity of the laser for two beams of different 
widths w (much thinner than the longitudinal size of the 
BEC cloud Wc)- The beam with w « Wc/50 (plotted with 
the symbol + in maroon) is able to extract more atoms 
than the beam with w ~ Wc/^b (plotted with the symbol 
X in blue). As it can be appreciated in the plot, the num- 
ber of atoms extracted decreases as the laser intensity is 
increased. This is due to the fact that the number of 
emitted solitons increases with the intensity of the laser, 
thus the number of atoms per soliton diminishes. 

In Fig. [5l we show three different numerical simula- 
tions employing beams of different intensities. Each fig- 
ure corresponds with the three points labeled in Fig. 2] as 
(a), (b) and (c). The vertical range corresponds to times 
from t = to t ^ For t = 1100/iy±_ the laser 

is set into motion along the z axis. The upper profiles 
show the condensate reservoir and the extracted atoms 
distributions. 



IV. ATOM EXTRACTION WITH MOVING 
TWEEZERS 

Another interesting possibility is to use a configura- 
tion of moving lasers with variable velocities. In this 
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FIG. 4: [Color online] Dependence of the percentage of ex- 
tracted atoms (N/Nbec) on the intensity of the laser beam 
power (P/Pmax) for two different widths of the outcoupling 
beam: w « Wc/50, plotted with the symbol + (in maroon) 
and w ~ Wc/25, plotted with the symbol x (in blue). The 
rest of the parameters are indicated in Fig. [5] 




FIG. 5: [Color online] (a-c) Examples of atom extraction 
from a condensate corresponding to different powers corre- 
sponding to the points labeled with (a), (b) and (c) in Fig|4] 
Shown are the density profiles for a time propagation of 
t = 1500/!/x. The laser is separated from its original posi- 
tion at t = llOO/:/^. (d-e) Pseudocolor plot indicating the 
full evolution of the BEC in the time span t G [0, 1500/!/±]. 



case, the extraction of atoms is obtained when the laser 
traverses the condensate, trapping particles through its 
path. The number of extracted atoms Ng varies with 
the velocity, and with the main parameters of the beam. 
We have analyzed, by means of numerical simulations, 
the dependence of on the beam intensity for several 
velocities in two different cases: the first one correspond- 
ing to a moving linear trapping potential and the second 
one corresponding to our nonlinear tweezers. Physically, 
the first situation consist of an optical dipole trap which 
does not change the value of the scattering length. We 
assume that in that case the laser which creates this lin- 
ear trap has the same width and depth as in the nonlinear 
case. The only difference is that nonlinear interactions 
are suppressed. Our purpose with this comparison is to 
evaluate the effect of nonlinear interactions in the extrac- 
tion procedure. To this aim, we have employed the same 



Schrodinger equation model as in the previous section. In 
all the simulations, a laser beam of width w is displaced 
from z^Otoz^Oata fixed given velocity u, extract- 
ing a fraction of atoms N^/Nb ec from the reservoir. The 
BEC reservoir is centered around z = 0. 
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FIG. 6: Top: comparison between the number of extracted 
atoms by nonlinear (continuous line) and linear (dashed line) 
tweezers for different values of the velocity. The y axis indi- 
cates the number of atoms trapped by the different optical 
tweezers normalized to the total number of atoms that form 
the initial BEC reservoir. The x axis represents the velocity 
of the tweezers when they traverse the condensate. Bottom: 
detail from the top plot for v = 8.0 mms~^ to v — 8.5 mms^^. 



1.0 
0.75 

0.5 
0.25 



5 



5 



0.0 

1.0 
0.75 

0.5 
0.25 

O.Q 





(a) 




















r 
















T 1 


1 

















1.5 



4.5 



w ()im) 















(b) 


















u 


1,1- 



1.: 



1.9 



2.1 



2.2 



V /V 

d' 



FIG. 7: Fraction of atoms extracted from the reservoir as a 
function of (a) the width w and (b) the intensity / of the 
optical tweezers. The velocity was fixed to « = 5 mm s~^ 
in all cases. The continuous lines represent the data for the 
nonlinear optical tweezers. The dotted lines correspond to 
the linear case. 

The simulations reveal that, in the nonlinear case, the 
fraction of extracted atoms Ns/Nbec depends crucially 
on the beam parameters, making the process highly con- 
trollable by changing the velocity, intensity or width of 
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the beam. In Fig. [6] we show the dependence of the per- 
centage of extracted atoms on the velocity of the beams. 
The black dashed line represents the values obtained with 
the linear tweezers. The continuous black line shows the 
dependence of the extracted atoms on the velocity in the 
nonlinear case. As it can be appreciated in the plot, the 
efficiency of the nonlinear optical tweezers is much higher 
than in the linear case. Another dramatic difference in 
both configurations is the presence of sharp variations 
of iVg at some velocities in the nonlinear configuration, 
allowing more control on the number of extracted atoms. 

In Figs. [7] we plot the variation of Ne/NsEC versus w 
[Fig. [Zlja)] and the depth of the dipole optical potential 
Vd, measured in units of Vb = hi>^/2 [Fig. Efb)]. In both 
cases, the data were obtained by fixing the velocity of the 
beams to u = 5 mm s~^. As in Fig. [5] the continuous 
lines refer to the nonlinear tweezers and the dashed lines 
to the linear ones. As it can be appreciated in the cap- 
tions, the behavior is similar to the variations observed 
in Fig. [SI The effect of nonlinear interactions is the ex- 
istence of sharp variations in the number of extracted 
atoms at certain values of the width and intensity of the 
laser tweezer. This adds extra control possibilities which 
are not accessible with linear traps. 

From the practical point of view the regions with com- 
plex variations of the number of particles as a function 
of the parameters (for example the range of potentials 
around Vd/Vo ~ 2.0 in Fig. [Zl^b)) are interesting in 
order to provide control on very few number of atoms. 
This effect opens new possibilities of achieving creation 
of macroscopic superposition of quantum states from ar- 
rays of Bose-Einstein condensates [solj. On the other 
hand, these regions are reminiscent of the fractal windows 
in chaotic scattering and resonances in soliton collisions 
typical of the interactions of nonlinear waves, which have 
received a lot of attention recently [4l| . 



V. CONCLUSIONS 

In summary, we have proposed a novel mechanism for 
extracting atoms from a BEG reservoir. Our system uses 
optical tuning of nonlinear interactions between atoms to 
extract them from the trap. 

By means of numerical simulations of the mean-field 
model equations, we have shown that this optically- 
induced spatial variation of the scattering length allows 
to control the number of atoms extracted and the po- 
sition of the outcoupled solitary wavepacket in a simple 
and robust way. 

We have also compared our nonlinear tweezer concept 
with the action of a linear potential as the one generated 
by ordinary laser beams in order to illustrate the crucial 
effect of nonlinear interactions in the process. We have 
also described chaotic scattering of solitons. 

Our results provide new ways to control the prepara- 
tion of BECs with a controllable number of atoms and 
are fully testable with current BEG experiments and can 
be easily generalized to systems with different atomic 
species. 
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